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ABSTRACT 

We investigate the relation between total X-ray emission from star-forming galax- 
ies and their star formation activity. Using nearby late-type galaxies and ULIRGs 
from Paper I and star-forming galaxies from Chandra Deep Fields, we construct 
a sample of 54 galaxies spanning the redshift range z w — 1.3 and the SFR 
range ~ 0.1 — lO'^Moyr"^. In agreement with previous results, we find that the 
Lx — SFR relation is consistent with a linear law both at z = and for the 
z = 0.1 — 1.3 CDF galaxies, within the statistical accuracy of ^ 0.1 in the slope 
of the Lx — SFR relation. For the total sample, we find a linear scaling relation 
Lx/SFR pa (3.5 ±0.4) x 1039(crg s-i)/(M0 yr"!), with a scatter of w 0.4 dcx. About 
~ 3/4 of the 0.5-8 keV luminosity generated per unit SFR is provided by HMXBs. We 
find no statistically significant trends in the Lx/SFR ratio with the redshift or star 
formation rate and constrain the amplitude of its variations by < 0.1 — 0.2 dex. These 
properties make X-ray observations a powerful tool to measure the star formation rate 
in normal star-forming galaxies that dominate the source counts at faint fluxes. 

Key words: X-rays: galaxies - galaxies: star formation - galaxies: starburst - X-rays: 
ISM - X-rays: binaries. 



1 INTRODUCTION 

High-mass X-ray binaries (HMXBs) and the hot ionized 
inter-stellar medium (ISM) are the main contributors to 
the total X-ray output of normal (i.e. not containing a lu- 
minous active galactic nucleus (AGN)) star-forming galax- 
ies. It is well established that the collective X-ray lumi- 
nosity of HMXBs well correlates with the star formation 
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The hot ionized ISM contributes about ^ 1/4 to the ob- 
served X-ray emission from late-type galaxies in the stan- 
dard X-ray band, its luminosity has been also shown to scale 
linearly with star formation rate (SFR) ((Grimes et al.ll2005l : 
iLi fc Wand |2Q12| : iMineo. Gilfanov. fc SunvaevI |2012| . here- 
after Paper II). Thus, it has been proposed that the total, 
integrated X-ray luminosity f rom star-forming galaxies can 
be us e d as a proxy of the SFR |Grimm. Gilfanov. fc SunvaevI 
l2003l : iRanalh. Comastri. fc Settil |2003| )." Although not en- 
tirely free from its own systematic uncertainties and con- 
taminations, the X-ray based SFR proxy is less affected 
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by the interstellar extinction and cosmological passband 
redshift, than conventional SFR indicators. Furthermore, 
the Lx — SFR scaling relation does not experience signif- 
icant cosmological evolution up to redshifts of z ~ 1 — 2. 
This has been initially suggested based on direct mea- 
surements of the Lx/ SFR r atios fo r several galaxies in 
Chandra deep fields ([Crimm. Gilfanov. fc SunvaevI 20031 : 
iLehmer et "al]l2010l : iMineo. Gilfanov. fc Sunvaevll2012l ) and 
was further supported by calculations of the maximal contri- 
bution of X-ray faint star-forming gala xies to the unresolved 
part of the Cosmic X-ray backg round (iDiikstra et al."20ia) 
and stacking analysis results (|Cowie. Barger. fc Hasingen 
I2OI2I ). These properties make the X-ray based SFR proxy a 
powerful tool to measure the star formation rate in distant 
galaxies. 

The most significant systematic effect which can com- 
promise the X-ray-based SFR measurements is the contam- 
ination by the emission of the central supermassive black 
hole, the AGN. Indeed, even low luminosity AGN, with 
log(Lx) ~ 42 can outshine a ~ lOOM© /yr starburst. As pop- 
ulations of bright gala xies are mainly composed of AGN (fo r 
recent results, see e.g. IXue et al. l l201ll : lLehmer et al.|[201^ ). 
SFR measurements using X-ray luminosity can be applied 
only to a relatively small fraction of bright galaxies. In this 
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case, a careful investigation of the nature of each galaxy 
is required in order to separate late-type from early-type 
galaxy populations. On the contrary, among faint sources, 
^ 10~^^ erg/s/cm^, the majority are star-forming galax- 
ies located at moderate and large redshifts, z ~ .5 — 3 
l|Cowie. Barger. fc Hasingedl2012| - TLehmer et al.ll2012l '). This 
makes the X-ray based SFR proxy a powerful tool to mea- 
sure the star formation rate in faint galaxies, where it can 
be u sed en masse, to infer the cosmic star formation history 
(e.g. ICowie. Barger. fc Hasingeill2012f ). 

The aim of this paper is to obtain the Lx — SFR scal- 
ing relation for the total X-ray luminosity and to investi- 
gate its behavior in a broad range of redshifts. For the red- 
shift z — 0, we use the sample of star-forming galaxies and 
ULIRGs (ultra-luminous infrared galaxies) from Papers I 
and II. We then select normal star-forming galaxies from 
the Chandra Deep Fields (CDF) expanding the local sample 
towards cosmologically interesting redshifts and high star 
formation rates. We combine these data in order to calibrate 
the Lx — SFR scaling relation in a broad range of redshifts 
and star formation rates. 

The structure of the paper is as follows. In Section 2 
we briefly summarize the selection criteria and properties of 
the local sample. In Section 3 we describe selection of late- 
type galaxies from the Chandra Deep Fields data and the 
procedures used to calculate their X-ray luminosities and 
star formation rates. The Lx — SFR relation is derived in 
the Section 4 and its redshift dependence is investigated in 
the Section 5. In the Section 6 we summarize our findings. 

Throughout this paper we assume fiat ACDM cosmol- 
ogy with Ho = 70 km/s/Mpc, = 0.3 and = 0.7. 
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Figure 1. The SFR — plane. Different types of galaxies are 
plotted with different symbols: Irregulars in (green) filled squares, 
spirals in (blue) filled triangles, ULIRGs in (red) empty circles, 
LIRGs in (orange) crosses, interacting systems in (grey) empty 
stars. Star-forming galaxies from Chandra Deep Fields -South 
and -North arc plotted as empty squares (purple) and filled cir- 
cles (gold) respectively. The dashed lines correspond to constant 
stellar-mass-to-SFR ratio. 



2 LOCAL GALAXIES 

2.1 The sample 

We constructed the local sample of star-forming galaxies by 
merging the sample of nearby galaxies resolved by Chandra 
from Paper II (Sect. 2, Table 1), and the sample of un- 
resolved LIRGs (luminous infrared galaxies) and ULIRGs 
defined in Paper I (Sect. 2.2, Table 2). The local sample 
contains 29 star-forming galaxies, their parameters are sum- 
marized in Table [T] 

The star formation rates and stellar masses (A/*) were 
determined in the Paper I, based on far-infrared, UV and K- 
band luminosities. These quantities were measured for the 
same spatial regions as used for computing the X-ray lumi- 
nosity, defined in the Sect. 5 of Paper I. The star forma- 
tion rates and stellar masses of the local sample span broad 
ranges, from ~ 0.1-290 Mq yr"^ and ~ 3x 10*-2x lO" Mq 
respectively. 

2.2 Total X-ray luminosity 

The data preparation was done following the standard 
CIACQ threads (CIAO version 3.4; CALDB version 3.4.1), 
exactly in the same way as described in Sect. 3 of Paper I. 

For resolved galaxies, listed in the first part of the Table 
[U the source spectra were extracted in the 0.5-8 keV energy 

^ http://cxc. harvard. edu/ciao3.4/index. html 



band and the background spectra were created as described 
in Sect. 3.3 of Paper II. Spatial regions for extracting source 
spectra are defined in the Paper I. Both source and back- 
ground spectra and the associated ARF and RMF files were 
produced by using specextract script and modeled using 
XSPEC V. 12.3. lx. The spectra were grouped in order to 
have minimum 15 counts per channel to apply the x fit- 
ting. For most of the spectra a good fit was obtained with 
a two component model: a thermal plasma (mekal) plus a 
power-law, corrected for the Galactic absorption. A second 
photoelectric absorption, accounting for the intrinsic absorp- 
tion and/or a second thermal component was required for 
some of the galaxies. It was decided based on the result of 
the F-test, using the probability of 10"'^ as the threshold (see 
Paper II for details). For each galaxy, we used the best fit 
model to calculate the count-to-erg conversion coefficient. 
Using this coefficient, the energy fiux was computed form 
the background subtracted source count rate. 

For the sample of unresolved starburst galaxies pre- 
sented in the second part of Table [1] we adopted 0.5-8 keV 
luminosities from the Paper I. 



3 STAR-FORMING GALAXIES FROM 
CHANDRA DEEP FIELDS 

In order to extend the Lx — SFR relation towards higher 
redshifts and star formation rates, we selected a sample of 
late-type galaxies in the Ch andra Deep Fields. W e started 
from the samples studied bv lLehmer et al.] (|2008l ') . who se- 
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Table 1. The local (z = 0) sample of star- forming galaxies. 



Galaxy 


D 


Ref._^ 


Hubble type 








log L*?*£ 




(Mpc) 








(10^° M,i 


) 


(ergs-1) 


RESOLVED GALAXIES 


NGC 0278 


11.812-4 


(2) 


SAB(rs)b 


4.1 


0.7 




39.52 


NGC 0520 


27.811? 


(2) 


Pec 


11.6 


4.7 




40.53 


NGC 1313 


4.1 ± 0.2 


(3) 


SB(s)d 


0.44 


0.1 




39.75 


NGC 1569 


1.9 ±0.2 


(4) 


IB 


7.8 X 10-2 


2.8 X 10" 


2 


38.17 


NGC 2139 


26.7l|j 


(5) 


SAB(rs)cd 


3.8 


0.91 




40.48 


NGC 3079 


18.2 


(1) 


SB(s)c 


6.0 


4.0 




40.30 


NGC 3310 


19.8 


(1) 


SAB(r)bc pec 


7.1 


0.98 




41.02 


NGC 3556 


io.7il:« 


(11) 


SB(s)cd 


3.1 


1.7 




39.87 


NGC 3631 


24.3 ± 1.6 


(7) 


SA(s)c 


4.6 


2.6 




40.79 


NGC 4038/39 


13.8 ± 1.7 


(8) 




5.4 


3.1 




40.57 


NGC 4194 


39 1+7!) 


(2) 


IBm pec 


16.8 


2.1 




40.92 


NGC 4214 


2.5 ± 0.3 


(9) 


IAB(s)m 


0.17 


3.3 X 10- 


2 


38.53 


NGC 4490 




(2) 


SB(s)d pec 


1.8 


0.39 




40.32 


NGC 4625 


8.2li;I 


(2) 


SAB(rs)m pec 


0.09 


7.0 X 10- 


2 


37.92 


NGC 5253 


4.1 ± 0.5 


(6) 


Im pec 


0.38 


5.9 X 10- 


2 


38.31 


NGC 5474 


6.8 


(10) 


SA(s)cd pec 


0.18 


9.1 X 10- 


2 


39.05 


NGC 5775 




(2) 


Sb(f) 


5.3 


6.3 




40.88 


NGC 7090 


7.6 


(1) 


SBc 


0.29 


0.22 




39.08 


NGC 7541 


34.9+^-? 


(5) 


SB(rs)bc pec 


14.7 


4.9 




40.39 


NGC 7793 




(11) 


SA(s)d 


0.29 


0.15 




38.53 


UGC 05720 




(2) 


Im pec 


1.8 


0.37 




39.77 



UNRESOLVED GALAXIES 



IRAS 17208-0014 


183.0 


(1) 


ULIRG 


289.9 


10.3 


41.40 


IRAS 20551-4250 


179.1 


(1) 


ULIRG 


139.4 


7.5 


41.63 


IRAS 23128-5919 


184.2 


(1) 


ULIRG 


139.6 


7.0 


41.88 


IRAS 10565-1-2448 


182.6 


(1) 


ULIRG 


156.8 


9.9 


41.42 


IRAS 13362-1-4831 


120.9 


(1) 


LIRG 


54.8 


17.3 


41.81 


IRAS 09320-1-6134 


164.3 


(1) 


ULIRG 


137.1 


13.1 


41.39 


IRAS 00344-3349 


84.0 


(1) 


LIRG 


20.1 


1.5 


41.23 


NGC 4676 


98.2 


(1) 




15.8 


6.6 


40.92 



" References for distances: (lllSanders et al.l ||2003|'), (2')ITu11vI |l988l '), (3')lGrise et all ll2008h . 
(4) iMakaroya fc Karachentsevl ||2003|'). (5) IWillick et all hgOTj). (61 ISaha et al. l ||2006| '). (71 
Ismail et al. 1 1|2005|1. (81 ISaviane. Hibbard fc RichI j2004ll. ( 9l lDrozdovskv et llT j2002ri . (10) 



Drozdovskvfcl<araclientse^ ~ il2000l 1 . (lU lTullv et al.l l l2(io9l 1. 



Star formation rate from Spitzer and GALEX data (see Sect. 6 Paper I). 
Stellar mass from 2MASS data (see Sect. 5 Paper I). 

Total X-ray luminosity of galaxies, i.e. gas and X-ray binary contribution, in 0.5—8 keV 
band, measured as described in Sect. [2^21 



lected late-type galaxies in the Chandra Deep Field-North 
(CDF-N) and Extended Chand ra Deep Field-South (E- 
CDF-S), and bv lXue et all (j201ll ). the 4 Ms Chandra Deep 
Field-South ( 4 Ms C DF-S) survey. For the South field, the 
iLehmer et al.l (|2008l ) sample was based on 1 Msec centered 
on the Chandra Deep Field-South and four flanking fields 
250 ksec each, while IXue et all |2oi3) covered only the 
Chandra Deep Field-South but to 4 times larger depth. 

These two samples were filtered according to the follow- 
ing criteria: (i) rigorous determination of the morphological 
type appropriate for a star-forming galaxy; (ii) availability 
of 0.5-8 keV X-ray flux measurement; (iii) availability of 
1.4 GHz radio flux mea surement for the SFR calculation or 
SFR determinations bv lLehmer et all (|2008l '): ( iv) the value 
of specific SFR (SFR/M*) larger than 10-^°yr-\ 



3.1 Selection by the morphological type 

The sample from iLehmer et al.l (|2008l ) includes 225 X-ray 
sources. Among them, 121 are classified as AGN candidates, 
while the remaining 104 sources are classified as normal late- 
type galaxies based on their rest-frame optical colors. For 
the aims of our analysis, we app lied two main selection cri- 
teria to the ILehmer et all (|2008l ') sampl e: we excluded i) all 
the AGN candidates and ii) objects detected only in one of 
the sub-bands, for which full band luminosity was not avail- 
able. After this selection, 52 objects were left, 37 from the 
CDF-N and 1 5 frorn the E-CDF-S. As it was pointed out by 
ILehmer et all (|2012l '). a single color division does not per- 
fectly disentangle the early-type and late-type galaxy pop- 
ulations. In order to exclude a possible early-type galaxy 
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Table 2. Late- type galaxies from the Chandra Deep Fields. 





Keasiiiit 


r lelQ 


Morpn. 
Type 


log Alif 

/AT \ 


iog -t^0.5-8kcV 
(erg s^i) 


'^1.4GHz 

( Nyj 


logivi.4GHz 
(erg s-i) 


or K 

(Mq yr ) 


J033206. 42-274728. 7 


1.02 


4Ms CDF-S 


3 


11.57 


41.80 


80.0 


30.55 


197.91 


J033211. 54-274713. 3 


0.58 


4Ms CDF-S 


6 


10.23 


41.08 


201.0 


30.35 


125.73 


J033217.90-275100.1 


0.12 


4Ms CDF-S 


4 


9.50 


39.52 


100.0 


28.55 


1.95 


J033218. 04-274718. 8 


0.73 


4Ms CDF-S 


8 


12.11 


41.59 


524.0 


30.92 


457.02 


J033220. 28-275222. 3 


0.34 


4Ms CDF-S 


6 


9.86 


40.69 


106.0 


29.60 


21.99 


J033228. 00-274639. 7 


0.25 


4Ms CDF-S 


6 


9.83 


40.01 


86.0 


29.16 


7.98 


J033229. 88-274424.5 


0.08 


E-CDF-S 


4 


10.05 


40.18 


1130.0 


29.24 


9.61 


J033229. 99-274404.8 


0.08 


E-CDF-S 


3 


9.88 


40.24 


482.0 


28.86 


4.07 


J033233. 00-275030. 2 


0.67 


E-CDF-S 


4 


10.64 


41.42 


42.0 


29.91 


45.52 


J033236. 18-274931.8 


0.55 


E-CDF-S 


8 


10.33 


41.53 


48.0 


29.65 


24.76 


J033237.27-275127.5 


0.69 


4Ms CDF-S 


6 


10.01 


41.23 


41.0 


29.80 


34.96 


J033238. 82-274956. 3 


0.24 


E-CDF-S 


3 


9.98 


40.48 


63.0 


28.99 


5.44 


J123559. 74-1-621550.3 


0.38 


CDF-N 


4 


10.56 


41.28 


212.0 


30.00 


55.25 


J123619. 47-1-621252.9 


0.47 


CDF-N 


3 


10.43 


41.31 


108.0 


29.92 


46.0 


J123634.47-I-621213.4 


0.46 


CDF-N 


8 


10.88 


41.53 


233.0 


30.22 


92.28 


J123632.45-I-621037.7 


0.14 


CDF-N 


6 


9.11 


39.08 






0.26^ 


J123634.53-I-621241.3 


1.22 


CDF-N 


6 


10.60 


42.39 


230.0 


31.20 


887.34 


J123646. 67-1-620833.6 


0.97 


CDF-N 


3 


10.89 


41.93 


80.1 


30.53 


188.27 


J123649. 72-1-621313.4 


0.47 


CDF-N 


6 


9.87 


41.14 


49.2 


29.56 


20.32 


J123651. 12-1-621031.2 


0.41 


CDF-N 


4 


10.23 


41.10 


95.0 


29.72 


28.89 


J123652. 83-1-621808.1 


0.25 


CDF-N 


3 


10.30 


40.53 


350.0 


29.80 


35.33 


J123653. 37-1-621140.0 


1.27 


CDF-N 


8 


10.44 


42.30 


65.7 


30.71 


285.27 


J123659. 92+621450.3 


0.76 


CDF-N 


8 


10.38 


41.45 


47.0 


30.05 


61.93 


J123708. 32-1-621056.4 


0.42 


CDF-N 


6 


10.07 


41.26 


45.1 


29.36 


12.64 


J123716. 82-1-621007.8 


0.41 


CDF-N 


3 


10.21 


40.89 


63.3 


29.55 


19.65 



R eferences fo r redsh ift and fields: E-CDF-S and CDF-N galaxies arc from lLehmer et"al ] l|200Sl '). 4Ms CDF-S galaxies are 
from lXue et~all ll201lh . 

Morphological Type from iBundv. Ellis, fc Conselicd ||2005| '). according to the following classification: 3=Sab, 4=S, 6=Irr, 
8=Merger. 

Determination from lLehmer et al. 1 1 I2OO8I) rescaled to Salpctcr IMF for 0.1 - 100 mass range. 



contamination of the lLehmer et all (j2008l l sample, we cross- 
matched the selected 52 sources with the morphological cat- 
aloRS of galaxies in the GOOD S North and South fields 
(|Bundv. Ellis, fc Conselic3l2005l ^. Using a match radius of 
1.5" we found a total of 44 pairs, 11 for the E-CDF-S and 
33 for the CDF-N sample. Two galaxies of the latter sample 
(J123603.26+621111.3 and J123627.32+621258.1) are classi- 
fied as early-type and were excluded. We also excluded eight 
galaxies without match in the GOODS catalogs, obtaining 
a final sample of 42 galaxies. 

As a check, we cross-correlated the obtained sample of 
CDF-N sources with the "High-SFR sample" in Paper I (Ta- 
ble 3), which includes seven galaxies from the Hubble Deep 
Field North (HDF-N). For seven HDF-N sources, we found 
six counterparts. The remaining galaxy, 123716.3+621512, 
was not present in the GOODS catalog. Correspondingly, 
we did not include it in our analysis. 

The 4 Ms C handra Deep Field-South (CDF-S) survey 
(|Xue et al.ll201ll ) includes 740 X-ray sources, of which 578 
are classified as AGN and 162 as normal galaxies. Follow- 
ing same strategy as above, we excluded all AGN candi- 
dates as well as sources for which broad band X-ray lu- 
minosity was not available. Late -type galaxies from the E- 
CDF-S, already included in the iLehmer et al.l (|2008l ') sam- 
ple, were not considered further. This shrank the sample 
to 94 galaxies. We cross-correlated these sources with the 



morphological catalog of galaxies in the GOODS South field 
iBundv Ellis, fc Conseli"c3l2005l ). Using a match radius of 
1.5" we found 75 matches. Among them we selected only 
sources classified as late-type or merging galaxies. Sources 
with no match were excluded. The obtained sample from the 
4 Ms Chandra Deep Field-South includes 48 star-forming 
galaxies. 

After this selection we were left with a sample of 90 
reliably identified star-forming galaxies from both Chandra 
Deep Fields. The sample will be reduced further after the 
cross-correlation with radio catalogs (see Sect. 



3.2 X-ray luminosity 

Both iLehmer et al] (|2008l ') and IXue et all (|201ll) provide 
rest-frame X-ray luminosities in the 0.5-8 keV band. In both 
cases the K-correction was p erformed assuming a power 
spectrum. ILehmer et al.l (|2008l ) used a photon in dex F = 2, 
appro priate for star-forming galaxies, whereas IXue et al.l 
((20 111 ) used a photon index of F — 1.8, more appro priate 
for AGNs. For sources from the ILehmer et all (|2008l ) sam- 
ple we used their original lumi nosity values. We recomputed 
luminosities for galaxies from IXue et al.1 (|201ll '). using their 
fiuxes and spectroscopic redshifts and adopting a photon 
index F = 2. 
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3.3 Star formation rate 

The infrared (IR) and ultraviolet (UV) luminosities are good 
proxies of the star formation activity in late-type galaxies. 
However, the use of these two estimators for CDF galaxies is 
limited by both the lack of rest-frame UV spectroscopic data 
and the sensiti vity limit of th e avail able IR observations. For 
these reasons, iLehmer et al.l (|2008l ) were able only to com- 
pute upper limits on SFR for most of the sources in their 
galaxy sample. To measure star formation rates in the CDF 
gala xies, we used 1.4 GHz radio emis sion. It has been shown 
(e.g. iBelll I2OO3I : ISchmitt et"aLl I2OO6I ') that this estimator is 
compatible with both IR- and UV-based estimators used 
for the local resolved galaxy sample (Sect. [2]). We cross- 
correlated the X-ray sample of star-forming galaxies pro- 
duced in the section mi with catalogs of radio sources based 
on VLA observ ations of CD Fs. For the CDF-N we us ed the 
iRichardj (|200(]| ') catalog, and lKellermann et all (|2008h cata- 
log for the CDF-S. Using a match radius of 1.5" we found a 
total of 24 matches with VLA sources: 12 from CDF-N and 5 
from E-CDF-S sampl es oflLehmer et al] (|2008l '). 7 from the 4 
Ms CDF-S sample of IXue et all (120111 ). The radio flux den- 
sities 51.4 GHz were converted to luminosities using the stan- 
dard procedure with K-correction. For the latter, we used 
radio spectral indices a as provided by VLA catalogues, ex- 
cept for 4 galaxies in CDF-N. For these sou rces we a s sume d 
a spectral index of a = 0.8, according to ICondonI (|l992h . 
T he radio-ba sed SFRs were estimated using the calibration 
of lBeiil (|2003l ): 

SFR(M0yr-') = 5.55 x 10-'''Li.4GHz(ergs-') (1) 

Among late-type galaxies without a radio counterpart 
in VLA catalogs, only one (J123632. 45+621037.7) has a n 
SFR determination in the catalog of iLehmer" et all (|2008h . 
For the rest of sources only upper limits on SFR were pro- 
vided. For the galaxy mentione d above we used the SFR es- 
timate bv ILehmer et al.) (|2008l ). which is based o n both IR 
and UV luminosities assuming a iKroupal (|200lD IMF. We 
rescaled the SFR to the Salpeter IMF (a factor of ~ 2.1) for 
0.1 — 100 mass range. 
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2. The - SFR relation. The (green) squares indi- 

cate galaxies from the local [z = 0) resolved sample, the (orange) 
triangles are LIRGs and ULIRGs from the unresolved sample, 
described in Sect. [2] and presented in Table [T] The filled (blue) 
circles are late-type galaxies from the Chandra Deep Field South 
and the and empty (red) circles are late-type galaxies from the 
Chandra Deep Field North, selected as described in Sect. [3] and 
presented in Table [2] The solid line shows the linear scaling re- 
lation, obtained using the whole sample of 54 galaxies and given 
by eq. 



catalog jGiavalisco et al.] 120041 ') . and z-band magnitudes 
fro m the GOODS-S Hubb le Space Telescope (HST) cata- 
log (jCiavalisco et al.ll2004l ). 

All the sources which reached this step satisfied to the 
SFR/A4 > 1 X 10"^'' yr"^ criterion. 



3.4 Stellar mass and specific star formation rate 

As a final check, we verified that selected CDF galaxies 
satisfy the selection criterion applied for the sample of re- 
solved galaxies in Paper I - that their specific SFR exceeds 
SFR/Af, > 10-^° yr-\ 

This was rat her straightforward for late-type galax- 
ies selected from ILehmer et all (120081 ) ■ for which stellar 
mass measureme nts were available. For CDF-S galaxies from 
IXue et al.l(|201ll ). we used the multi-wavelength optical data 
provided in the same catalog. We calculated the stellar 
masses from the K-band luminosity and (r — z) color calibra- 
tion o f the mass-to- light ratio from lZibetti. Chariot, fc Rixl 
(|2009l ): 

In particular, we derived Lk/Lk,q based on K-band 
AB magnitudes from the MUSYC catalog (|Tavlor et al.l 
[iooi), and computed (r — z) colors based on AB mag- 
nitudes in R-band from the ESQ 2.2 m WFI R-band 



3.5 The final CDF sample 

The final CDF sample consists of 25 galaxies from the Chan- 
dra Deep Fields, 13 objects from CDF-North and 12 from 
CDF-South. They are listed in Table [21 along with their 
parameters. 

The entire sample, including 21 local (resolved) galax- 
ies, 8 LIRGs and ULIRGs and 25 CDF galaxies (54 galaxies 
in total) is shown in the SFR-M* plane in the Figjl] 



4 ix - SFR RELATION FOR TOTAL 
LUMINOSITY 

We approximate the Lx ' — SFR data (in Fig|2]) with the log- 
Unear model log Lx ' = log if + /3 log SFR. This model was 
applied to local sample and CDF sample separately. Both 
slope P and normalization K were set as free parameters of 
the fit. To fit the data we used least square minimization. 
The results are given in Table O which demonstrates that 
the best fit parameters for the local and CDF samples are 
compatible. We therefore performed the joint fit to the whole 
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Table 3. Summary of the parameters for Lx — SFR relations obtained from 
least-squares fit. 





FREE SLOPE 




LINEAR FIT 


Sample 


log/\: 


/3 


cr 


logK 


cr 








(dex) 




(dox) 


2 = galaxies 


39.53 ± 0.09 


1.04 ±0.08 


0.42 


39.56 ±0.08 


0.42 


CDF galaxies 


39.62 ± 0.15 


0.95 ±0.09 


0.33 


39.54 ±0.07 


0.33 


all 


39.55 ± 0.08 


1.00 ±0.05 


0.38 


39.55 ±0.05 


0.38 



Note. The parameters are relative to the least-squares fit to the data with the 
relation log Lx = log K + j3 log SFR, respectively setting the slope /3 free and 
fixing it to unity, a is the dispersion around the best-fitting relation. See Sect. 
12.21 and 13.21 for the definition of the X-ray luminosities, Sect. 6 of Paper I and 
l3.3l of present work for the definition of SFR for z = sample and CDF galaxies 
respectively. 



sample of star-forming galaxies that includes both local and 
CDF sources, using same model. We obtained a slope /3 = 
1.00 ± 0.05, suggesting a linear relation. The linear scaling 
relation obtained by fixing the slope to unity is shown in 
Fig. [2] and it is given by the following equation: 

^o°5-8kcv(ergs'') ^ (3.5 ±0.4) x lo'' SFR (Mq yr"^) (3) 

where the uncertainty was computed from the scatter of 
the points around the best fit. The results of all fits are 
summarized in Table [3] 

The scale factor in eq. @ is larger than the scale factor 
obtained for the collective luminosity of HMXBs in resolved 
galaxies (eq.(22) Paper I) by a factor of ~ 1.34. Therefore, 
we can conclude, that the unresolved emission contributes, 
in average, ~ 1/4 of the total X-ray luminosity of star- 
forming galaxies generated per unit of SFR. This result is 
consistent with the preliminary estimate obtained in the Pa- 
per I based on the comparison of the Lx — SFR relations of 
the HMXBs in the resolved galaxies and total emission from 
unresolved galaxies. It is also consistent with the scaling re- 
lation for the ISM emission obtained in the Paper II. 

In the conclusion of this section we note that the re- 
sult of the fitting of the Lx — SFR relation depends on 
whether the fit was performed in the logarithmic or lin- 
ear space. In the linear case, it also depends on whether 
the least square fit in the form Lx = A x SFR was per- 
formed or individual Lx/SFR ratios were averaged (note 
that in the log case it is same). Note also that large scatter 
of points, much bigger than statistical uncertainties, ren- 
ders the minimization less appropriate. Doing the least 
square fit to the linear quantities, we obtain the scale of 
Lx/SFR f» (2.7±0.2)xl0^^ (erg s-^)/(M0 yr"^). Averaging 
the Lx/SFR ratio in the linear space, we obtain yet another 
number: Lx/SFR « (5.0 ± 0.6) x 10^^ (erg s-^)/(M0 yr'^), 
i.e. nearly twice larger than in the eq.([3|. The value ob- 
tained in the logarithmic fit is smaller, because it gives equal 
weights to the low and high points0 As this seems to be a 
more correct approach and in order to be consistent with 
the Paper I, we used the fit in the log space as default. 



^ Indeed, for two measurements yielding yi = a and y2 = 10 a 
the logarithmic (geometrical) average is {y)igg ~ 3.16 a while the 
linear average is (i;);^^ = 5.5 a, i.e. is nearly twice larger. 



4.1 Comparison with previous results 

iRanalli. Comastri. fc Settil (|2003l 'l studied the Lx - SFR re- 
lation using X-ray data from ASCA and BeppoSAX satel- 
lites. In the Paper I, we converted their X-ray luminosity 
and SFR to make them consistent with the analogous quan- 
tities used here and derived: Lo.s-skoV ~ 1.6 • 10=*^ X SFR, 
i.e. more than ~ 2 times smaller than our eq.Q. The 
rather large discrepancy is likely caused by the confusion 
in the definition of SFR proxies (see detailed discussion in 
the Paper I). On the other hand , the Lx — Lib. relation of 
iRanalli. Comastri. fc Settil (|2003l ). transformed to be com- 
patible with our definitions is Lo.s-skeV ~ 2.4 • 10~* x Lir. 
The scale in this relation is « 27% larger that the scale 
in the corresponding relation for HMXBs from the Paper I 
(eq.(23)), consistent with ~ 1/4 contribution of hot ISM and 
other unresolved emission components to the total luminos- 
ity - see the discussion after the eq.Q above. This proves 
that the underlying Lx — I/ir relations are fully compatible 
and the discrepancy in the Lx — SFR relation is due the 
difference in the use and definitions of SFR proxies. 

IPersic fc Rephaelil |20o3) investigated the relation be- 
tween the total X-ray emission from normal and starburst 
(ULIRGs) galaxies in the 2-10 keV band and the SFR de- 
rived from infrared luminosity. Using their eq. (10) and con- 
verting the 2-10 keV luminosity to the 0.5-8 keV band (a 
factor of 1.28) we obtain: Lo.s-gkov ~ 4.9-10^^ x SFR. This 
relation is consistent, within a factor of ~ 1.4 with our eq. 
((3]). On the other hand, in the Paper I we found a rather 
large disagreement between our Lx — SFR relations for the 
HMXB luminosity. 

iLehmer et al] (|2010l ') recently studied the relation be- 
tween 2-8 keV luminosity and SFR for a sample of nearby 
LIRGs using Chandra observations. They obtained the hard 
band luminosities by modeling the X-ray spectra in 0.5- 
8 keV band using Galactic absorption, a thermal compo- 
nent and a power-law component. This is consistent with 
the model we used to fit the X-ray spectra of our sample of 
LIRGs and ULIRGs and determine their X-ray luminosities. 
In order to c ompare our results w ith the Lx — SFR relation 
obtained by iLehmer et al] (|2010l ). we converted their 2-8 
keV luminosity to the 0.5-8 keV band. We assumed same 
model as above, with average quantities from our ULIRGs 
spectral fits (F ~ 1.7, kT ~ 0.6 keV, A^h ~ 2 x 10^^ cm"^ 



© 2010 RAS, MNRAS 000, [HH] 



Total X-ray emission from star-forming galaxies 7 



,40.5 



40 



0) 



"39.5 



CO 



O 



39 



38.5 



"I 1 1 1 1 1 1 \ 1 r 

. - CDF-S 
~. - CDF-N 



■ 

▲ 


■ 

▲ 


-t — ' 






0) 



,40.5 



40 



'39.5 



O 



39 



38.5 



TTTTj 1 — I I I I I ll| 1 — I I I I I ll| 1 — I I I I lllj 1 — I I I I III 

. - local (z = 0) 
~ _ CDF-S ~ 
. . - CDF-N 




iiilJ I I I I I I I 



0.1 1 10 100 1000 

SFR [Mg yr-i] 



Figure 3. Dependence of the L^?' — SFR relation on redshift (left panel) and star formation rate (right panel). The symbols without 
error bar indicate individual galaxies. The (black) filled circles with error bars show their average values within each redshift (left) or 
SFR (right) bin. The vertical error bars are uncertainties of the average values computed from the rms of points in the given bin. In 
the left panel, the horizontal dot-dashed line shows the average local L^*/SFR ratio. Its uncertainty is shown by the narrower (orange) 
shaded area, while the wider (grey) shaded area shows the rms of the points. On the right hand panel, the horizontal line and the shaded 
areas show similar quantities computed using all galaxies. 



and thermal-to-power-law normalization ratio of ~ 0.5). For 
the linear (fir s t line in their Table 4) Lx — SFR relation in 
iLehmer et all (|201(]| ') we obtain a scale factor Lx/SFR ~ 
2.9 X 10^® (ergs-^)/(M0yr-^). This result is also in a rea- 
sonable agreement with our eq. 



5 CONSTRAINTS ON THE REDSHIFT AND 
SFR DEPENDENCE OF THE Lx SFR 
RELATION 

The composition of the sample (29 galaxies at the redshift 
2 = and 25 galaxies at z ~ 0.1 — L3) allows us to constrain 
a possible redshift evolution of the X-ray emission produced 
per unit SFR. To this end, we grouped the CDF galaxies 
into three redshift bins of the equal width spanning from 
z = ^ \,o z = 1.3 and computed the average Lx*/SFR ratio 
and its uncertainty within each bin. Same quantities were 
computed for the local sample of galaxies, to characterize 
the 2 = population. The result of this calculation is shown 
in Fig. |3l It is evident from the figure, that the average 
values for CDF galaxies are consistent, within uncertainties, 
between each other and with the average quantity for local 
galaxies. To confirm this, we applied the Spearman's rank 
test and we found no statistically significant correlations, 
neither for the individual values not for the average points, 
with the corresponding probabilities of P --^ 27% and P ^ 
20% respectively. 

Thus, the total X-ray emission produced per unit SFR 
by late-type galaxies does not show statistically significant 
trends with the redshift up to z ~ 1.3. The sensitivity of 



this analysis is sufficient to eliminate the changes larger 
than ^ 0.1 — 0.2 dex (by factor of ^ 1.3 — 1.6), although 
as a caveat we should note that the high redshift bin con- 
tained only 4 galaxies. This conclusion is in agreement with 
results of the X-ray stacking analysis of late-type galaxies 
iLehmer et al.|[2008h and with calculations of the maximal 
contributi on of star-forming ga laxies to the unresolved CXB 
intensity (|Diikstra et al.ll20ld ). 



Recent population synthesis calculations by 
iFragos et al.l |2012i ) predict that the ratio of the 0.5- 
8 keV luminosity to the SFR increases by a factor of 
~ 1.3 — 1.4 between redshift 2 = and z ~ 1.3. Such an 
evolution is within the uncertainties of our results. 



Finally, in the right hand panel of the Fig[3]we constrain 
the possible dependence of the Lx/SFR ratio on the SFR. 
Using both local and CDF galaxies, we bin the data into 
3 SFR bins and compute the average and its uncertainty 
(rms-based) for each redshift bin. The formal ^ = 4.8 for 
the binned data, with 2 dof. The probability of having such 
(or bigger) value of the ^ is ~ 9%, which is equivalent to 
~ 1.5 a detection. This may or may not indicate the exis- 
tence of some SFR dependence in the Lx — SFR ratio. If 
any, the amplitude of this dependence also does not exceed 
~ 0.15 dex (a factor of ~ 1.4). Finally, we note that due 
to the ob server's bias, the effect of statistic s of the small 
numbers jGilfanov. Grimm, fc Sunvaevll2004l ) is not strong 
in our sample. 
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6 SUMMARY 

Based on the sample of nearby resolved galaxies more distant 
ULIRGs at intermediate distances and star-forming galax- 
ies from the Chandra Deep Fields, we construct a sample 
of 54 star-forming galaxies spanning the range of redshifts 
from z = up to z = 1.3 and the range of star formation 
rates SFR ~ 0.1 - 10^ Mo/yr (Fig[T]|. Using this sample, 
we calibrate the Lx — SFR relation for the 0.5-8 keV band 
luminosity (Fig[5]|. We find that Lx — SFR dependences for 
the local and CDF samples are consistent with linear rela- 
tions with the typical accuracy of ~ 0.1 in the slope. The 
linear Lx — SFR relation obtained for the entire sample is 
given by the eq.Q. We did not find any statistically signifi- 
cant trends in the scaling relation with the redshift and star 
formation rate with the upper limit on the possible varia- 
tions in the Lx/SFR ratio of ~ 0.1 — 0.2 dex (a factor of 
~ 1.3 — 2.6) (Fig[3|. This property makes the X-ray emission 
a powerful tool to measure star formation rate in a broad 
range of redshifts and star formation regimes, which can be 
applied en masse to faint distant galaxies. 
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